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ABSTRACT 






This thesis presents a comparison of three Magnetic 
Anomaly Detection (MAD) models: a c ro s s - c or r e 1 a t i on 

detection model, a square law detection model, and a model 
referred to as the OPTEVFOR detection model. FORTRAN and 
BASIC programs for the three detection models are included 
in this thesis. The programs yield detection probabilities 
for straight line encounters. Magnetic signal values for 
the straight line encounters are an additional output. 
Plots of lateral range curves and magnetic signal values are 
presented. A discussion of the required parameters is 
included in the thesis to facilitate the use of the 
programs. The parameters that were considered in the 
comparison of the three detection models are: magnetic 

noise, aircraft and submarine headings, submarine 
displacement, and the vertical separation between submarine 
and aircraft. 
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I. 



INTRODUCTION 



This thesis presents a comparison of three Magnetic 
Anomaly Detection (MAD) models. The comparison is in terms 
of probabilities of detection that were computed using the 
models. Two of the models, the cross-correlation model and 
the square law model, have been used to model sonar 
detection [Ref. 1: pp. 343-357]. The third model, referred 
to as the OPTEVFOR model, is a slant range threshold 
detection model. The results of the comparisons are 
presented in graphical and tabular form. In addition, plots 
of magnetic signals for selected lateral ranges and noise 
levels are shown. The effects of noise, aircraft and 
submarine headings, submarine displacement, and vertical 
separation are also indicated. 

The models were implemented using the FORTRAN and BASIC 
programs^ that are listed in Appendix A. For those 
interested in using the programs for other investigations, 
an input parameter discussion is provided in Chapter 3. To 
use the FORTRAN program, the user specifies the input 
parameters in an input file. After execution of the 
program, an output file is generated that contains 

1 The programs are based on an unpublished BASIC program 
by R.N. Forrest for an H.P.- 85 microcomputer. 
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probabilities of detection for each of the three models. In 
addition, magnetic signal values and magnetic signal values 
plus random magnetic noise values for one of the encounters 
generated by the program are included in the output file. 
An IBM GRAFSTAT graphical package was used to produce the 
graphics in this thesis. 

To use the BASIC program, the user must interactively 
enter the input parameter values for each encounter. After 
execution of the program, an optional hardcopy printout 
supplies the input parameter values and a table of detection 
probabilities for each of the three models (see Appendix A). 
Following this, lateral range curves are displayed to the 
user for immediate observation. A typical program run 
producing 21 detection probabilities for each model requires 
approximately 10 minutes of computing time on an Atari 800 
microcomputer . 



II. 



MODEL DESCRIPTIONS 



A. THE CROSS-CORRELATION AND SQUARE LAW DETECTION MODELS 

The cross-correlation and square law detection models 
are described in detail by Forrest [Ref. 2: pp. 33-35 ]. The 

models are based on the following assumptions: the noise is 

gaussian, and the signal sample points are such that 
adjacent magnetic noise samples are independent. 

The magnetic signal values, as measured by an 
aircraft's magnetometer for the cross -cor re la ti on and 
square law detection models, are the submarine magnetic 
field values in the direction of the earth's magnetic field 
at the positions of the magnetometer. The submarine's field 
is assumed to be a dipole field, and the aircraft and the 
submarine are assumed to keep constant speeds and headings 
during an encounter. 

For the cr os s- correl at ion model, a complete prior 
knowledge of the magnetic signal is required. 
Operationally, this suggests that a signal file, which 
contains a replica of the signal for each possible encounter 
situation, would be required. The model describes a 
perfect detection system with respect to the noise model 
that is used. For the square law model, a signal replica is 
not required. This model might be considered to describe 
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the limiting detection capability for an automatic system 
that does not use information about the shape of the 
magnetic signal. 



B. OPTEVFOR MAD DETECTION MODEL 

The OPTEVFOR model is described by Forrest [Ref. 3: pp . 
7-8]. In characterizing the submarine magnetic signal as a 
simple dipole signal, the U.S. National Defense Research 
Committee, [Ref. 4: p. 20], reports that the magnetic signal 
of the submarine "varies as the inverse cube of the distance 
from the source". In an OPTEVFOR report [Ref. 5: p. 1, end. 
1], the results of a regression analysis on empirical peak 
to peak signal output against slant range between submarines 
and aircraft are reported. These results also suggested 
this inverse cube relationship for the magnetic signal. 
This relationship is the basis for the OPTEVFOR detection 
model . 

The model has a deterministic mode and a stochastic 
mode, each of which involves the following parameters: the 

submarine magnetic moment (M), an Operator Recognition 
Factor (ORF), the average peak to peak magnetic noise (N) in 
the operating area, and a slant range (R). The relationship 
between these quantities is given by: 

c M 

(ORF) N 



R= 



1/3 



( eqn 2.1) 
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The value of the constant c is 0.10 for M in oersted 
centimeters^, R in meters, and N in gamma. 

In the deterministic mode, detection occurs if and only 
if the aircraft's slant range from the submarine at CPA is 
less than or equal to R. This mode yields a rectangular 
("cookie cutter") lateral range curve with the probability 
of detection equal to 1 for an encounter where the slant 
range at CPA is less than or equal to R, and 0 when it is 
greater than R. 

The stochastic mode allows a more uncertain approach to 
detection by allowing a gradual rise in probability of 
detection as the slant range at CPA decreases. In this mode 
one sets the probability of detection at R equal to 50 
percent, and the lateral range curve is given by p^ = (J)(x); 
where it is understood that ^ is the standard normal 
cumulative distribution function and x is determined by the 
following equation: 

R - CPA 
x = 

(AL) R . ( eqn 2.2) 

In this equation, CPA is the magnitude of the slant range 
distance at CPA, and R is the calculated range from Equation 
2.1. The product (AL)R represents a standard deviation. 
The value of AL can be considered to be determined by "the 
combined uncertainty and variability in the values of M, N, 
and ORF" [Ref. 3: p. 8]. Two values of AL (.20 and .01) are 
shown in Figure 2.1. If empirical data was available, the 

1 3 



P(DH) 



value of AL could be chosen to provide a best fit to the 
observed results. Note, as AL approaches 0, the stochastic 
mode approximates the deterministic mode. 



al - 0.20 





-1000 0 1000 
METERS 



Figure 2.1. 



Lateral Range Curves for Different Values 
of AL . 
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III. 



INPUT PARAMETERS 



The input parameters for the FORTRAN program are all 
contained in one input file. This allows parameter values 
to be easily changed without recompiling the main program or 
subroutines. Also, with a few changes, this program could 
be altered to operate in conjunction with a larger program 
to yield a probability of detection on an individual MAD 
run. 

The input parameters are divided into four areas for 
discussion. They are: (1) sample interval, (2) earth 

magnetic field, (3) submarine moments, and (4) other 
inputs . 

A. SAMPLE INTERVAL 

The choice of a sample interval is discussed by Forrest 
[Ref. 2: pp. 27-30]. In the program, the total observation 
time in seconds over which the samples are taken is entered 
in T7. This time should be long enough to encompass a 
"complete signal" at the maximum expected detection slant 
range . 

As the slant range from the submarine to the 
magnetometer increases, the distance over which a 
significant magnetic signal is present at the magnetometer 
also increases. Figure 3.1 graphically shows the difference 



in the amount of time that a signal is present for slant 
ranges of 200 meters and 805 meters. In this thesis, the 
total time for a straight line encounter is assumed to be 20 
seconds. As can be be seen from Figure 3.1, a 20 second 
interval adequately covers the significant portion of the 
magnetic signal for an 805 mmeter slant range at CPA. 



200 METERS CPA SLANT RANGE 




805 METERS CPA SLANT RANCE 




Figure 3.1. Magnetic Signals for Slant Ranges of 200 Meters 
and 805 Meters. 

The time between samples is set equal to the reciprocal 
of twice the upper bandpass filter frequency of the MAD 
sensor. A value of 0.9 Hz was suggested for use by Texas 
Instruments [Ref. 6: p. 112] as an upper bandpass filter 
limit in a discussion on the effects of noise on a MAD 
system. This value yields a time interval between samples 
of 0.55 seconds. 

The sample interval length and the false alarm rate (the 
expected number of false alarms per hour) determine the 
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false alarm probability. The false alarm rate (F2) is 
assigned a value of 3 based on a report by OPTEVFOR [Ref. 5: 

p. 2.1]. 

B. EARTH MAGNETIC FIELD 

Input values for the earth magnetic field intensity and 
inclination, or dip angle, may be taken from two Defense 
Mapping Agency Hydrographic Center charts, [Refs. 7 and 8 
respectively], or approximated by using a program. If chart 
values are entered, the earth field intensity must be in 
units of gamma and the inclination in decimal degrees. The 
program used to determine the intensity of the earth field 
and inclination is based on a simple dipole field model that 
is described by Forrest [Ref. 9: pp. 39-43]. 

Table III-1 displays the program output values of 
inclination in decimal degrees and earth magnetic field in 
gamma for selected geographic locations. In addition, 
corresponding values obtained from the Defense Mapping 
Agency Hydrographic Center Charts Number 30 and Number 39 
are also displayed. The last three columns are the average 
slant range in meters at which a 50 percent probability of 
detection is obtained for the three program detection 
models. The program input parameters for these slant ranges 
were the same as the base case, except for the following 
differences: a sample interval time of 40 seconds, aircraft 
and submarine headings of 0 degrees, and a submarine 
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one oersted = 10 D gamma 



displacement of 7,000 tons. The correlation between the 
slant ranges, comparing the chart values and model values, 
was found to be 95 to 96 percent for the three models. This 
suggests that, even though differences exist between the 
chart values and model values, there is a high degree of 
correlation in the final output. 

A limitation to the simple dipole field model is that 
it does not give an angle of declination (variation) with 
sufficient accuracy. As a result, all headings entered 
into this program must be in magnetic degrees. The Phoenix 
Corporation [Ref. 10: pp. 24-25] reports on geomagnetic 

field models that can represent the earth field "with 
overall accuracies better than approximately 150-200 gammas 

o 

in magnitude and .2 in direction of the field." This 
degree of accuracy is not needed for this program, but a 
simplified version of one of these models that provided 
satisfactory angles of declination would be beneficial if 
the program were to be incorporated into a larger model that 
utilized true headings as inputs. 

C. SUBMARINE MAGNETIC DIPOLE MOMENT 

If a submarine's magnetic dipole moment is known for the 
geographical location and the submarine's magnetic heading, 

^Private communication from R.N. Forrest, who 
investigated the use of the simple dipole model for this 
purpose . 
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the following values may be entered in the program: (1) P, 
its magnitude in oersted centimeters cubed, (2) A, its 
direction in decimal degrees relative to magnetic north, and 
(3) B, its depression angle from the horizontal in decimal 
degrees. If it is not known, these values must be 
calculated for a specific location and magnetic heading. A 
program is included in the main program that can be used to 
calculate these values. The program is based on a model 
described by Forrest [Ref. 9: pp. 35-38 ]. The input to the 
program is submarine displacement in tons. The program also 
contains coefficients which relate displacement to magnetic 
moment. The values used in the program are based on values 
cited by Texas Instruments [Ref. 6: p. 4], 

The past history of the submarine is represented by the 
permanent longitudinal, transverse, and vertical moments of 
the submarine (M4, M5 and M6 in oersted centimeters cubed). 
For the examples in this thesis, it was assumed that 
effective deperming had been performed and program default 
values of zero were used. 

D. OTHER PARAMETERS 

1 . Headings and Speeds 

Since the simple dipole earth field model used by 
the program does not produce accurate angles of declination, 
magnetic headings are required. In addition, the headings 
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must be in decimal degrees. The input parameters for 
submarine speed and aircraft speed are entered in knots. 

2 . Noise 

The magnetic noise is assumed to be such that 

adjacent magnetic noise samples are independent. This 

assumption is based in part on the filtering that is 

performed on the magnetic signal by the processing system in 

a MAD detection sensor. The standard deviation of the noise 

in units of gamma is the value entered into SI. This value 

can be approximated from operational data by taking from 

one-fourth to one-sixth of the measured peak to peak 
I 

magnetic noise. [Ref. 2: pp. 28-29] 

The OPTEVFOR detection model incorporates a value of 
average peak to peak magnetic noise (N) in the inverse cube 
law calculation. In the program, the value of N is 
determined by multiplying the SI entry by four. 

3 . Operator Recognition Factor ( ORF ) 

The ORF is the value of the ratio of magnetic signal 
to magnetic noise for which the average operator would 
detect a signal 50% of the time in the presence of 
background noise for a false alarm rate of 3 per hour. An 
ORF value of 3 was suggested for use by OPTEVFOR [Ref. 5: p. 

4.12] . 

4. Distance Parameters 

Two parameters, R8 and N7, are used to define the 
points plotted on the lateral range curves. R8 is the 
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maximum positive value of the lateral range in meters for 
which a lateral range curve value is to be computed. N7 
represents the number of lateral range curve values that are 
to be computed from the maximum lateral range to zero 
lateral range. 

The vertical separation (Z) is the sum of the submarine 
depth and aircraft altitude in meters. 
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IV. RESULTS 



Program outputs of the three models for a set of base 
case conditions are presented in this section. Outputs for 
variations from the base case are also presented. The 
lateral range of an encounter (the horizontal separation 
between the submarine and magnetometer when the magnetometer 
is at CPA) for a 50% probability of detection is used as a 
measure of comparison. Signal and signal plus "noise" 
traces for several cases are presented. The traces are 
based on the signal and noise models that are part of the 
cross-correlation and square law models. These idealized 
signal traces appear to have the characteristics of actual 
signal traces. This suggests that the signal and noise 
models might be used for training purposes. 

A. BASE CASE 

The base case conditions are listed in Table IV-1. The 
table is ordered in the same manner that the values are read 
into the program. An annotation of each entry is included 
for clarity. 

Figure 4-1 presents the lateral range curves for the 
base case. Points on the lateral range curves are indicated 
by the first letter of the name of the model from which they 
were derived. The slight asymmetry of the cross-correlation 
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detection model and square law detection model curves is 
reflective of the shape of the signals that are 'received' 
in these models. 



Table IV- 1 . Input Parameters for the Base Case 



1 .8 
20.0 

3.0 
0 

30.0 

60.0 

45.0 

10.0 

315.0 

220.0 
0 

0 

0 

4000.0 

200.0 
0.1 

1500.0 

50.0 

3.0 
0.2 
0025 



twice the upper bandpass limit in seconds 
sampling time interval in seconds 
false alarms per hour 
Enter inclination (1 = yes, 0 = no)? 
area of operation latitude in decimal 
degrees 

area of operation longitude in decimal 
degrees 

submarine magnetic heading in decimal 
degrees 

submarine speed in knots 

aircraft magnetic heading in decimal 

degrees 

aircraft speed in knots 

Enter submarine moment (1 = yes, 0 = no)? 
Enter earth field (1 = yes, 0 = no)? 

Enter submarine perm moments (1= yes, 

0 = no ) ? 

submarine displacement 
vertical separation in meters 
noise (standard deviation) in gamma 
maximum lateral range in meters 
divisions of lateral range 
ORF (Operator Recognition Factor) 
variability factor for OPTEVFOR model 
lateral range iteration number for the 
magnetic signal and signal plus noise in 
the output file 



Table IV-2 lists lateral detection ranges and 
corresponding slant detection ranges at CPA for a 
probability of detection equal to 50 percent for the cross- 
correlation and square law detection models. An equivalent 
ORF value for each model is also listed. Due to the 
asymmetry of the lateral range curves for the cross- 
correlation and square law models, the average of the two 50 
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Figure 4-1. Lateral Range Curves of the Cross-Correlation 
(C), Square Law (S), and OPTEVFOR (0) Models for the Base 
Case . 



Table IV-2. The Lateral Detection Ranges, Slant Detection 
Ranges, and ORF's of the Three Models for the Base Case. 





Lateral Detection 
Range (meters) 


Slant Detection 
Range (meters) 


ORF 


Cross -Cor relation 


885 


907 


.21 


Square Law 


685 


71 4 


.44 


OPTEVFOR 


31 8 


376 


3 



percent detection ranges was used as the lateral detection 
range. The equivalent ORF values for the cross-correlation 
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and square law detection models were calculated using the 
slant detection range values with the following equation, 
which was obtained from Equation 2.1: 

c M 

ORF = 

N . eqn 4.1 

For the base case, the magnitude of the submarine field (M) 
at the submarine is 6.35 x 10° orested- cm , the noise (N) is 
.4 gamma, and the value of the constant (c) is .1. This 
suggests that, in order to detect a magnetic signal 50 
percent of the time with a false alarm rate of 3 per hour, 
the magnetic signal to magnetic noise ratio should be .21 
for an ideal cross-correlation detector and .44 for an ideal 
square law detector. 

Using the ORF values, the cross -corre la ti on and square 
law detection models can be used to describe the performance 
of an operator. To do this, a modified value of the 
standard deviation (o') of the input noise can be used. The 
modified value is equal to (ORF ) (o') / .2 1 for the cross- 
correlation detection model and ( ORF) (o')/. 44 for the square 
law detection model. With these modifications, the two 
models can be used to describe the detection capability of 
an operator with a specified ORF. An example of a lateral 
range curve with the modified noise standard deviation for 
an ORF of 3 is presented in Figure 4.2 for each model. 
These curves are comparable to the lateral range curve for 
the OPTEVFOR model that is shown in Figure 4.1. 
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CROSS-CORRELATION: NOISE - 1.43 GAMMA 



SQUARE LAW; NOISE - 0.68 GAMMA 



D 

o 

ST 





Figure 4.2. Cross-Correlation and Square Law Lateral Range 
Curves to Describe the Performance of an Operator with an 
ORF of 3. 



The 


automatic 


MAD system 


manu 


f a ctured 


by 


Canada 's 


CAE 


Electronics Ltd. 


is expected 


to p 


roduce a 


50 


pe rcent 


in - 


cr ea se 


in detect 


ion slant range 


[Ref. 11 ]. 


Using 


the 



detection slant range for the OPTEVFOR model of 376 meters, 
a 50 percent improvement would yield a detection slant range 
of 564 meters. The ORF for a detection system with this 
capability would be .88. The cross-correlation and the 
square law detection models could be used to yield lateral 
range curves for a system with an ORF of .88 by using a 
noise standard deviation equal to .88 (O')/. 21 and .88 ( 0-* ) / 
.44 respectively. Figure 4.3 shows the lateral range curves 
of the two detection models with a 50 percent improvement in 
slant range detection. Note, with the modified noise 
standard deviations, the models are essentially equivalent 
for the cases considered. 
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CROSS— CORRELATION; NOISE » .2 GAMMA 



SQUARE LAW; NOISE - 2 GAMMA 



13 



Q 

BT 





Figure 4.3. The Cross -Correlation and Square Law Models to 
Describe LRC ' s for the CAE Automatic Detection System. 



SIGNAL AT CPA OF 0 METERS 



SIGNAL + NOISE AT CPA OF 0 METERS 





Figure 4.4. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise at a Lateral Range at CPA of 0 Meters for the 
Base Case. 



Figures 4.4 and 4.5 present the magnetic signal and a 
representation of magnetic signal plus magnetic noise that 
would be received under the base conditions by a 
magnetometer with a lateral range of 0 meters and of 780 
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meters. The signal plus noise trace was generated from 
signal plus noise values obtained by adding a signal value 
to a gaussian noise value. The gaussian noise value was 
generated by multiplying the standard deviation of the input 
noise by a pseudo normal random number from a population 
with mean 0 and variance 1. The pseudo normal random numbers 
were generated using LLRANDOM II, a resident program at the 
Naval Postgraduate School computer [Ref. 12: p. 2.2]. 

SIGNAL AT CPA OF 780 METERS SIGNAL + NOISE AT CPA OF 780 METERS 





Figure 4.5. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise at a Lateral Range at CPA of 780 Meters for 
the Base Case. 

The magnitude of the magnetic signal shown in Figure 4.4 
is very large in comparison to the background noise. The 
peak to peak signal to noise ratio is approximately 14 to 1. 
An operator would have little difficulty identifying the 
signal in this signal plus noise trace. 

Conversely, the magnetic signal shown in Figure 4.5 is 
small compared to the background noise. The peak to peak 
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signal to noise ratio is .35. The probabilities of 
detection for the lateral range of 780 meters are: .95 for 
the cros s- corr el at ion detection model, .28 for the square 
law detection model, and 0 for the OPTEVFOR detection model. 
It seems apparent that an operator would have a difficult, 
if not impossible, time in detecting this signal at a 
reasonable false alarm rate. 

B. DIFFERENT NOISE INPUTS 

The first variation on the base case shows the effect of 
different noise inputs. The standard deviation ( O' ) of the 
peak to peak noise is the input parameter that is varied. 
Table IV-3 lists the different O’ values and the 
corresponding lateral detection ranges. 

Table IV-3. The Effect of Noise on Detection Range. 

Standard Deviation Lateral Detection Range in Meters 

of Noise in 

Gamma Cross- 





Correlation 


Square Law 


OPTEVFOR 


.005 


2 2 5 0 ( 2259) * 


1 7 9 2 ( 1 803) * 1 000 


( 1 020) 


.01 


1832 


(1843) 


1 446 


(1460) 


782 


(807) 


.05 


1110 


(1128) 


868 


(890) 


427 


( 472) 


.1 


885 


( 907) 


685 


(714) 


31 8 


(375) 


.5 


51 2 


( 550) 


382 


(431 ) 


90 


(219) 


*The 


numbers in 


parentheses 


are the 


slant ran 



distances in meters. The vertical separation is 200 meters. 

Figure 4.6 displays lateral range curves for the three 
models when the standard deviation of the noise is .01 
gamma. These three curves show an increase in lateral 
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detection range over the base case. Note that the asymmetry 
of the cross-correlation and square law detection models is 
more apparent in Figure 4.6 than it was in Figure 4.1. 

Figure 4.7 displays the magnetic signal (which is the 
same as the signal in Figure 4.5) and the magnetic signal 
plus magnetic noise at a horizontal distance of 780 meters 
when the magnetometer is at CPA. The signal to noise ratio 
is 3.5. The figure suggests that a MAD operator, in this 
case, should have the ability to detect a signal at 780 
meters lateral range with a satisfactory false alarm rate. 




Figure 4.6. Lateral Range Curves for the Three Models with 
the Standard Deviation of the Noise Set to .01 Gamma. 
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SIGNAL AT CPA OF 780 MEIERS 



SIGNAL + NOISE AT CPA OF 780 METERS 





Figure 4.7. Magnetic Signal and Magnetic Signal Plus 
Magnetic Noise with the Standard Deviation of Noise = .01 
Gamma at 780 Meters Lateral Range. 

C. DIFFERENT HEADINGS 

The headings of a submarine and an aircraft in an 
encounter have an effect on detection ranges. The effect of 
different headings was investigated using the square law 
detection model, and the results in terms of lateral 
detection ranges are presented in Table IV-4. This table 
suggests that a submarine should choose a magnetic heading 
of either East or West, and, for an encounter, an aircraft 
should also choose a magnetic heading of East or West. 

Figure 4.8 shows lateral range curves for a submarine 
heading North and an aircraft heading East. In this case, 
both the cro ss -cor re la ti on and square law detection model 
lateral range curves display noticeable asymmetry. The 
OPTEVFOR detection model lateral range curve is symmetric 
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Table IV-4. Square Law Lateral Detection Ranges for 
Different Submarine and Aircraft Magnetic Headings 



Aircraft 
Headings 
( magnetic ) 





0 


45 


0 


700 


650 


45 


724 


637 


90 


730 


654 


1 35 


712 


685 


180 


685 


640 


225 


71 2 


624 


270 


730 


646 


315 


724 


685 



Submarine 



90 


135 


180 


498 


640 


686 


505 


670 


71 3 


527 


654 


730 


505 


636 


724 


498 


650 


700 


501 


672 


724 


51 9 


646 


730 


501 


624 


71 3 



Headings (magnetic) 



225 


270 


31 5 


640 


498 


650 


624 


501 


676 


646 


524 


646 


682 


502 


624 


650 


498 


640 


636 


506 


672 


654 


533 


654 


682 


505 


636 




Figure 4.8. Lateral Range Curves for the Submarine Heading 
North and the Aircraft Heading East. 



33 



but, like the curves for the other models, it shows an 
increase in detection ranges over those for the base case 
(where the submarine is heading NE and the aircraft is 
heading NW ) . 

The APAIR MOD 2.6 [Ref. 13: p. 83] simulation uses a MAD 
detection model that accounts for the change in a 
submarine's magnetic moment (which is dependent on changes 
in submarine heading) by using a parameter labeled DFACTR 
(degradation factor for heading). In the model, D (a 
modified slant range at CPA) determines the probability of 
detection. The value of D is determined using the following 
relation : 

D = DC ( 1 - DFACTR x A), eqn . 4.1 

where DC is the slant range at CPA and A is the acute angle 
in decimal degrees between the submarine heading and an 
East-West bearing. The probability of MAD detection is 
determined from a table of probability of detection against 
slant range. A uniform (0, 1) random number is drawn to 

determine whether or not the submarine is detected. The 
average slant detection ranges (computed from Table IV-4, 
where the vertical separation is 200 meters) for submarine 
headings of North and East are 741 meters and 545 meters 
respectively. These ranges yield a value of .003 for 
DFACTR. The average slant detection range from Table IV-4 
for a submarine heading of NE is 682 meters; however, the 
slant range determined by a modified slant range of 545 
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meters and a DFACTR = .003 is 643 meters. If sin A instead 
of A is used in Equation 4.1, then DFACTR is .26 5 and the 
slant detection for a submarine heading NE is 670 meters. 
Since this is only a single data point and there is no 
supporting operational data, the modification is not 
proposed as one that should be adopted. However, this 
cursory analysis does indicate a way in which the programs 
presented in this thesis might be used by others. 

Table IV-5 lists lateral ranges for P(det) equal to 50 
percent for 3 submarine/aircraft heading combinations. The 
cross-correlation and OPTEVFOR detection model results show 
the same relationship as the results of the square law 
detection model. 

Table IV-5. Lateral Ranges for P(det) = .50 in Meters for 
the Three Detection Models. 



Submarine 

Aircraft 


45 

315 


0 

90 


90 

0 


Cross -Cor relation 


885 


934 


754 


Square Law 


685 


730 


498 


OPTEVFOR 


318 


358 


230 



For the detection ranges reported by OPTEVFOR [Ref. 5: 
p. 5.1], the effect of different headings was averaged out. 
That is, measurements were taken from the 16 possible 
combinations of the 4 cardinal submarine and aircraft 
headings in equal numbers and then averaged to yield an 
average slant detection range. But, as shown in Tables IV-4 
and IV-5, the models show significant variability in lateral 
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detection range for different submarine and/or aircraft 
headings . 

Figure 4.9 is included to show the lateral range curves 
when the submarine is headed East and the aircraft is headed 
North. These lateral range curves give the minimum lateral 
detection ranges for the different heading combinations. 
Also, for the cross-correlation and square law detection 
models, the lateral range curves are fairly symmetric. 




Figure 4.9. Lateral Range Curves for the Submarine Headed 
East and the Aircraft Headed North. 
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D. SUBMARINE DISPLACEMENT 

The submarine magnetic dipole moment program within the 
main program is used to calculate a submarine's induced 
magnetic moments. The program is based on a model described 
by Forrest [Ref. 9: pp. 35-38]. The model requires 

submarine displacement as an input. Table IV-6 displays 
results when the submarine displacement is doubled in each 
succeeding entry. 



Table IV-6. Slant Detection Ranges in Meters for 
Different Submarine Tonnages. 



Displacement Signal 


Slant Detection Ranges 


in Meters 


in tons 


Magnitude 
in oersted 
cm J 


Cross- 

Correlation 


Square 

Law 


OPTEVFOR 


1 000 


1 . 59x1 0 8 


590 


463 


236 


2000 


3.1 7x1 0 8 


732 


575 


297 


4000 


6.35x1 0 8 


907 


71 4 


376 


8000 


1 .27x1 0 9 


1 1 27 


885 


472 


1 6000 


2.54x1 0 9 


1 402 


1099 


597 


32000 


5.08x1 0 9 


1 724 


1363 


753 


As can 


be seen from 


column two in 


Table IV-6, the dipole 


moment is 


proportional 


to the displacement. Since 


the three 


detect ion 


models give a slant detection range 


that is 



proportional to the cube root of the dipole moment, doubling 
the submarine displacement should multiply the slant 

1 / O 

detection range by 2 ' (1.26). This is confirmed by 

comparing the slant detection ranges between the entries in 
Table IV-6. Doubling the displacement multiplies the slant 
detection range by 1.24 for the cross-correlation and square 
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law detection models and, as expected, by 1.26 for the 
OPTEVFOR detection model. 



Table IV-7 lists the displacement in tons of selected 
Soviet submarines. The values were taken from Combat Fleets 
of the World 1 982 /1 983 [Ref. 14: pp. 602 -61 4 ]. This table 
Table IV-7. Selected Soviet Submarine Displacements. 
Class Displacement in Tons 



Typhoon 
Delta III 
Yankee 
Echo II 
Victor I 
Charlie I 
Tango 
Foxtrot 
Whiskey 



25-30,000 

10,500-13,250 

8.000- 9,600 

5.000- 6,000 
4,300-5,1 00 

4.000- 4,900 

3.000- 3,700 
1 ,950-2,400 
1 ,080-1 ,450 



is presented solely for the purpose of the information it 
contains. The submarine magnetic dipole moment program 
should not be expected to give accurate estimates of these 
submarine's induced magnetic moments, since the program uses 
a value that relates displacement to magnetic moment that is 
based on submarines of World War II. 



E. VERTICAL SEPARATION 

Figure 4.10 shows three lateral range curves for a 
vertical separation of 500 meters. The OPTEVFOR detection 
model lateral range curve shows only a slight detection 
probability even when the aircraft passes directly over the 
submarine. The cross-correlation and square law detection 
model lateral range curves show an increase in lateral 
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detection range over the base case. The dip in the lateral 
range curves, for each of these models, suggests the complex 
variation of the magnetic signal with lateral range. 




Figure 4.10. Lateral Range Curves for a Vertical Separation 
of 500 Meters. 

Table IV~8 lists the lateral detection ranges for 
different vertical separations. It should be kept in mind 
that these values are for a single geographic location; 
consequently, they may not be representative of other 
locations. Mote that both the cross-correlation and square 
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law detection models lateral detection ranges increase with 
an increase in vertical separation until about 500 meters. 



Table IV-8. Lateral Detection Ranges for Different Vertical 
Separations . 



Vertical Separation 


Lateral Detection Range 


in Meters 


in meters 


Cross- 


Square Law 


OPTEVFOR 




Correlation 






1 00 


804 


61 4 


360 


200 


885 


685 


31 8 


300 


942 


720 


225 


400 


974 


724 


-- 


500 


980 


699 




600 


974 


629 




700 


936 


262 




* 

No longer attains a probability of 


detect ion 


equal to 50 


percent . 








A factor related to 


vertical separation is the effect of 


ocean wave noise on a 


MAD system. 


As the altitude of a 



magnetometer is decreased, the magnitude of the ocean wave 
noise increases. Because of the rate of this increase, for 
a given submarine and submarine depth there is a minimum 
altitude at which an aircraft should prosecute a submarine 
using MAD. Further investigation using an ocean wave noise 
model might be valuable. 
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V. CONCLUSIONS 



This thesis has presented a comparison of three MAD 
detection models. The cross -cor re la ti on detection model, 
which models an optimum detector under the conditions of the 
detection model, yields the maximum detection range for a 
set of given conditions. The square law detection model 
does not describe an optimum detector under the conditions 
of the model and yields shorter detection ranges. In the 
stochastic mode, with an appropriate choice for the 
parameter AL that determines the standard deviation, the 
lateral range curves for the OPTEVFOR detection model become 
similar to the other two detection models. Detection ranges 
for the OPTEVFOR detection model depend on the choice for 
the Operator Recognition Factor (ORF). With a value of 3 
for the ORF, it yields the shortest detection ranges. 
Adjusting the magnetic noise level by an amount proportional 
to the effective ORF, the cross-correlation and square law 
models can be used to describe the performance of an 
operator or an automatic detection system. 

The magnetic signal and magnetic signal plus noise 
traces appear to have the characteristics of actual signal 
traces. This suggests that the signal and noise models, 
which are the basis for the cross-correlation and square law 
detection models, might be useful for training purposes. 
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Variations on a set of base case parameters were used to 
show relative changes in the detection models. The 
parameters included: magnetic noise, submarine and aircraft 
magnetic headings, submarine displacement, and vertical 
separation. Significant results were the large asymmetry of 
the lateral range curves under certain conditions and the 
variation of the magnetic signal as shown by the changes in 
vertical separation. 

The FORTRAN and BASIC programs, along with an input 
parameter discussion, are included to facilitate the use of 
the three MAD detection models as they are implemented by 
the programs. 
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10 DIM G (200) » D1 < 100) » D2( 100) » K( 100) » X0( 100) 

IS DIM D5C100) 

20 DEG 

30 PRINT "MAX FREQ" 7 
40 INPUT FI 

S0 PRINT : PRINT "MAX FREQ = "?F1 
60 Tl-l/Fl 

70 PRINT "INTERVAL TIME =* " 7 : INPUT T7 
30 G-T7/2/T1 
30 H®INT(G) 

100 H»H+INT(2*<G-H)> 

110 M«2*H+1 

120 IF M> 200 THEN 70 

130 PRINT "INT TIME - " 7T7 

140 T7=T1*M 

150 PRINT "ADJ INT TIME = " 7T7 
1S0 PRINT "SAMPLE SIZE » "7M 
170 PRINT SPRINT "F/A RATE "7 
130 INPUT F2 

130 PRINT : PRINT "F/A RATE * " 7F2 
200 P 1 =F2* < M- 1 ) *T 1 /3630 
210 PRINT "PF a " 7 PI 

220 PRINT SPRINT "INPUT DIP ANGLE ( 1=YES» 0=NO) " 7 : INPUT A 
230 IF A=0 THEN 260 

240 PRINT "DIP ANGLE PHI " * s INPUT F 

230 GOTO 420 

260 DEG sLl=7S:L2=100 

270 PRINT SPRINT "LATITUDE "7 

230 INPUT L 

230 PRINT ".PRINT "LONGITUDE "7 
300 INPUT 0 

310 PRINT SPRINT "LAT =» " 7L s PRINT "LON =* "70 
320 F*SIN<0-L2)*C0S<L> S G=CQS(G-L2) *COS CL) :H=3INCL) 

330 U=G : V=H 
340 GOSUB 1300 

350 J=J-(30-Ll) :G=K*SIN(J) :H=K*COS< J) 

3E0 U=G:V=F 



370 GOSUB 1300 

330 F=K s G=0 s R=s- (COS <L1)*SIN(J) ) :Q=- CCQSCLl > *CQSC J) ) 

330 U=H:V=F 
400 GOSUB 1300 

410 F=ATN C2*( SIN ( J) /CCS( J) ) ) 

420 PRINT "PHI = " 7F 

430 PRINT SPRINT "DIPOLE COURSE " 7 s INPUT ClsPRINT "DIPOLE SPEED " 7 s INPUT VI 

440 PRINT SPRINT "SENSOR COURSE "7 s INPUT C2: PRINT "SENSOR SPEED " 7 s INPUT V2 

450 PRINT "DIPOLE COURSE - "7C1 SPRINT "DIPOLE SPEED « "7V1 

480 PRINT "SENSOR COURSE * " 7C2s PRINT "SENSOR SPEED =» " ?V2 

470 Wl— V2*SIN (C2) — V1*SIN (CD :W2=V2*C0S (C2) -V1*C0S (Cl ) 

480 U*W1 :V=U2 
430 GOSUB 1300 
500 C0=JSW0=K 



510 D3=UI0*T1*4. S3/3 

520 PRINT "REL COURSE => "7C0:PRINT "REL SPEED = " 7UI0 

530 PRINT SPRINT "INPUT DIPOLE MOMENT Cl -YES? 0=NO) " 7 s INPUT AA 

540 IF AA=0 THEN 5S0 

550 PRINT SPRINT "MAGNITUDE P " 7 s INPUT P: PRINT "HOR ANGLE w "7 
5E3 PRINT SPRINT "VERT ANGLE OMEGA " 7 s INPUT 3 
573 GOTO 340 

583 PRINT SPRINT "INPUT EARTH FIELD ( l«YESf 0-NQ!) INPUT AA 



INPUT A 
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533 IF AA=3 THEN 620 

633 PRINT : PRINT " EARTH FIELD ,# 7 : INPUT El 
613 GOTO 633 

623 E 1 =73000/ SQR < 3*C0S < F> *COS < F) + 1 ) 

633 PRINT SPRINT "EARTH FIELD = "?E1 
643 ri4=0:n5=3:n6=0 

653 PRINT SPRINT "INPUT PERM MOMENTS <l=YES,0=NO) INPUT AA 
6E3 IF AA=0 THEN 690 

673 PRINT SPRINT "LONG MOMENT "? s INPUT M4s PRINT "TRAN MOMENT INPUT MS 
6S3 PRINT "VERT MOMENT "7 : INPUT MS 

653 PRINT SPRINT "LONG MOMENT = "7M4: PRINT "TRAN MOMENT = "SMS: PRINT "VERT MOMEN 
T =■ "SMS 

700 Kl-7.35K2=1.6sK3*1.6 

713 PRINT SPRINT "DISPLACEMENT "S 

723 INPUT N1 

733 PRINT SPRINT "DISPLACEMENT = "SN1 
733 NN1=N1 

74Q MS=6 1 *K3*N 1 *S I N < F > s M3=M9+MS 

733 MS=E1^C0S(.F^>^N1^<K1>hCQS<:C1)‘*‘CCSCC1)+K2>»*SINc:C1)^SIN<C1)) 

763 Ml=M4*SIN<Cl)+M5*C0S<C:) :M2=M4*C0S<C1 )-M5*SIN<Cl) 

773 M7=El**<.C0St;F^i:Kl-K2^Nl'HSIN<Cl)^C0S(:Cl) ) 

7S3 Ml=M7+Ml:M2=MQ+M2 
793 U=Ml : V=M2 
200 GQSUB 1230 
213 A=J:lJ=M3:V=K 
223 GOSUB 1203 
233 P*K : B=J 

243 PRINT SPRINT "P = "SPsPRINT "w = "SAsPRINT "OMEGA = "SB 



233 X=P1 



263 GOSUB 1223 
273 V6=*Y 



220 V7=M*< 1— 2/9/M+Y*SQR < 2/9/M) ) A 3 



220 PRINT SPRINT 
900 PRINT SPRINT 
910 PRINT SPRINT 
920 INPUT SI 
925 S2=S1*4 
930 PRINT SPRINT 
243 PRINT SPRINT 
933 PRINT SPRINT 
263 PRINT SPRINT 
964 PRINT SPRINT 
2ES PRINT SPRINT 



"VERT SEPARATION "s: INPUT 2 
"VERT SEPARATION = "72 
"NOISE "S 



"NOISE = " SSI 

"MAX LATERAL RANGE " S : INPUT RS 
"NUMBER QF INCREMENTS "s: INPUT N7 

"MAX LATERAL RANGE = "SRS SPRINT "NUMBER OF INCREMENTS 
" ORF " S s INPUT ORF 
" ORF =■ " SQRF 



973 D4=R3/N7sNS=2*N7 

974 PRINT SPRINT " ALPHA " S s INPUT AL 
976 PRINT SPRINT " ALPHA =» "SAL 

923 L9=-RS 
993 FOR E=0 TO N8 
1300 X0=L2:X0CE)=L3 
1010 GOSUB 1620 
1020 L3=L3+D4 
JL030 GOSUB 1200 
1040 NEXT E 
1044 GOTO 1290 
1030 GRAPHICS 3s COLOR 1 
1060 XX=* INT <310/N8) 

1065 X0(0)=0 

1070 FOR 1=0 TO N8 

1080 di<:i)=intc<i-oi<:i))**i60> 

1090 D2Cn = INTCCl-D2<I))^l60) 

1095 D5CI)=INTCC1-Q5<I))*163) 

1130 x0a+n=x0d)+xx 

1110 NEXT I 

1120 PLOT X0(0).D1(0) 

1130 FOR 1=1 TO NS 
1140 DRAWTO X0CI)tD1CI) 



"7N7 
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1123 

1160 

1170 

11S0 

1120 

1132 

1134 

1135 
1138 
1200 
1220 
1220 
1235 
1330 
1310 
1323 
1323 
1335 
1340 
1333 
1360 
1370 
1330 
1320 
1400 
1413 
1423 
1430 
1440 
1450 
1460 
1470 
1420 
1430 
1300 
1505 
1513 
1515 
1323 
1525 
1533 
1535 
1343 
1545 
1550 
1560 
1605 
1610 
1620 
1S30 
1640 

1642 

1643 
1645 



NEXT I 

PLOT X0C0),D2C0> 

FOR TO N3 

DRPUTQ X0a),D2CI) 

NEXT I 

PLOT X0<0>, D5<0) 

FOR 1=1 TO N8 
DRAUTQ X0CI)fD5(I) 

NEXT I 

PRINT "PD FOR X FROM "7-RST" TQ " ? RS 
GOTO 1610 

PRINT “FOR HARD COPY ENTER ’ 1’ " ? : INPUT CC 

IF CCOl THEN GOTO 1050 

LPRINT '‘MAX FREQ = "?F1 

LPRINT “ADJ INT TINE = " T T7 

LPRINT "SAMPLE SIZE = " T M 

LPRINT "F/A RATE = "TF2 

LPRINT "PF = "TP1 

LPRINT "LPT = "1L 

LPRINT "LON = "TO 

LPRINT "PHI = " *F 

LPRINT "DIPOLE COURSE =» "TCI 

LPRINT "DIPOLE SPEED = "TV1 

LPRINT "SENSOR COURSE =» " T C2 

LPRINT "SENSOR SPEED = " T V2 

LPRINT "REL COURSE = " T C3 

LPRINT "REL SPEED = "TW0 

LPRINT "EARTH FIELD = "tEI 

LPRINT "LONG MOMENT = "Tf14 

LPRINT "TRAN MOMENT = " T MS 

LPRINT "VERT MOMENT = " T MS 

LPRINT "DISPLACEMENT = " TNN1 

LPRINT "P = "TP 

LPRINT "w = "TP 

LPRINT "OMEGA = " » B 

LPRINT "VERT SEPARATION = " T Z 

LPRINT "NOISE = "TS1 

LPRINT "MAX LATERAL RANGE = “ »RS 

LPRINT "NUMBER OF INCREMENTS * "TN7 



LPRINT : LPRINT "LTR RNG PDCCC) 

L2=»- R8 

FOR 1=0 TQ N8 

LPRINT L3» " "TD1 (I) *" " » D2< I) T " 

LS=LS- fc *D4 

NEXT I 

GOTO 1053 

PRINT "END" 

END 

U=X0: V=Z 
GOSUB 1300 
D=J:H0=K 

RH= <!0. 1*P/C0RF*S2) ) *0.333 

SIG=AL*RH 

X=< RH—H0) /SIG 



PD <SL) 



PD <QPT) ' 



"TD5CI) 



1646 GG SUB 1390 

1647 D5 (E) =Y 

1630 30=COSCB)*COSCC3-^) : J0=*ZQS<D)^COSCB)^SINCC3-A)-SIN<D)^SINCB) 

1660 N 0 =-<SIN<D)> 4 «CaSCB)^SINCC 0 -A) J-COSCD) *SINCB) 

1670 31=COSCF)^CaS<:C0) :Jl=»CQSCD)^CaS<F)^SINCC3)-SIN<D)^SINCF) 

16S0 Nl=-<SIN<D)^CaSCF)’^SINCC3))-CaSCD)>*‘SINCF) 

1630 K1=P/10/H0a3 

1700 A2=2*80*B1-J0*J1-N0*N1 : A1=3*<:n0*01+00*N1) : A0=2*N0*N1-B0*B1- J0*J1 
1710 S0=0 

1723 FOR 1=0 TO M-l 
1730 S 3 < I- <M—1 ) /2) *Q3 
1740 Q=S/H0 
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1753 G=1 / < 1+Q*Q) A 2. 5 
1763 G=< A2* ! Q*Q+A1*Q+A0) *G 
1773 G < I ) =G : S0=S0+G *6 
1730 NEXT I 
1730 RETURN 
1300 K2=SQR<S0> 

1310 K<E>=K2 

1320 VS=- V6+K 1 *SQR ( S0 ) / S 1 

1330 L0=K1*K1*S0/ <S1*S1> : A3=M+L0:B3=1+L0/ CM+L0) 

1343 V9=-SQ R < 2*V7 / B3 ) +SQR < 2* A3/ B3- 1 ) 

1350 X=V3 : GQSUB 1990 
1360 D1CE)=Y 
1370 X=V9 : GOSUB 1390 
1330 D2CE)=Y 
1390 RETURN 

1300 K=SOR<U*U+V*V> : IF K=0 THEN J=0:RETURN 
1905 UK=U/K:VK=V/K 

1907 IF UK) 0.999999 AND VK) 0.999999 THEN J=0: RETURN 

1903 IF UK) 0.999999 THEN J=-ATN VK/SQR C-VK*VK+1 ) )+90 : RETURN 

1909 IF VK) 0.999999 THEN J=0: RETURN 

1910 f1M=ATN <!UK/SQR<!— UK*UK+1 ) ) : J=— ATN< VK/SQR<!— VK*VK+1 } )+90 5 IF MM<0 THEN J=360-J 
1915 RETURN 

1920 Y=X : IF X)0. 5 THEN Y=l-Y 
1930 Y=SQR<LOG Cl/Y/Y)) 

1340 G0-2. 515517 : Gl=0. 302353 : G2=0. 010323 

1950 HI = 1- 432733 :H2=0. 139269 :H3=1 . 30SE-03 

1 960 Y=Y- < G0+ Y* < G 1 +G2* Y ) ) / < 1 +Y* C H 1 +Y* c H2+H3*Y > ) ) 

1970 IF X) 0. 5 THEN Y—Y 
1980 RETURN 

1330 Y=X:IF X <0 THEN Y*-Y 
2000 U»1/C1+0.2316413*Y> 

2010 Q1=0. 31933153:02=^0. 356563732:03=1. 73147793 : Q4=-l . 32122559 : Q5=l - 33027442 
2020 IF Y) 24. 23 THEN Y=0:GGTO 2070 
2025 PI =3. 14159265 

2030 Y=EXP<-<Y*Y/2>)/SQR<2*PI>*W*<Ql+W*<Q2+W>KQ3+W*<Q4+W*Q5) ) ) ) 

2070 IF X>0 THEN Y=l-Y 
2075 Y=< I NT C 10000*Y) ) / 10000 
2030 RETURN 



55 



MAX FREQ - 1.3 

ADJ INT TIME = 20.55555555 

SAMPLE SIZE = 37 

F/A RATE = 3 

PF = 0. 0166666666 

LAT = 30 

LON = 60 

PHI = 59.40076973 
DIPOLE COURSE = 45 
DIPOLE SPEED = 10 
SENSOR COURSE = 315 
SENSOR SPEED = 220 
REL COURSE = 312.397439 
REL SPEED =220.227153 
EARTH FIELD = 52506.5513 
LONG MOMENT = 0 
TRAN MOMENT - 0 
VERT MOMENT = 0 
DISPLACEMENT = 4000 
P = 634694392 
w = 32.63750751 
OMEGA = 27. 11173179 
VERT SEPARATION = 200 
NOISE = 0. 1 

MAX LATERAL RANGE = 1500 
NUMBER OF INCREMENTS = 15 



LTR RNG 


PDCCC) 


PD <SL) 


PD (.OPT) 
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